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Abstract 

O 
O 

CN | Angular momentum is important concept in physics, and its phase space 

properties are important in various applications. In this work phase space 
■ analysis of the angular momentum is made from its classical definition, and 

by imposing uncertainty principle its quantum properties are obtained. It is 
shown that kinetic energy operator is derived, but it has different interpreta- 
tion of its parts than in the standard treatment. Rigid rotor is discussed and 

i — i 

' it is shown what is its phase space representation. True rigid rotor is defined 

and also its phase space properties are discussed. 

ON 
O 

I. INTRODUCTION 

Angular momentum is a very important concept in dynamics of particles and within 
quantum mechanics its properties are very well understood |1J. One example where the 
theory has direct application is the rigid rotor model, which is basic for understanding rota- 
tional spectroscopy and collisions of molecules. These processes are described by quantum 
dynamics, but there are circumstances when classical dynamics is used as alternative. For 
example, rotational cross sections for two colliding molecules in principle can be calculated 
from quantum mechanics, but often it is a challenging task. Classical mechanics, on the other 
hand, is relatively simple to use, but there several problems in its implementation, say to 
calculate atom-molecule collision cross sections. The basic one is how the initial conditions 
are selected and the final results analyzed. For example, if one says that a molecule is in 
the rotational state with the quantum numbers I = 3 and m = 2 then the question is what 
to choose for the initial orientation and angular velocity which adequately represents it? 
Analogous problem had been analyzed for collisions involving only vibrational energy ex- 
change in atom-molecule collisions, where it was shown how to select initial conditions that 
adequately represent a particular vibrational state of the molecule. A recipe was suggested 
from which transition probabilities were successfully calculated || from classical mechanics, 
and by that it is meant that initial conditions were selected from a prescribed phase space 
density but dynamics is calculated from classical equations of motion. 

The problem is therefore how does one chooses a set of initial positions and velocities 
for a particle so that on average its angular momentum has given I and m values? In other 
words, what is needed is a function fi, m (r,p) with a property that the following averages are 
obtained 
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< (fx p) 2 >= J d 3 r d 3 p (rxp) 2 f l)m (r,p) = h 2 l(l + 1) (1) 

< fx p >= / d 3 r d 3 p fxp fi tm (f, p) = hm z 



In classical mechanics the function fi, m {f,p) is interpreted as a statistical weight from which 
the position and momentum are selected. It should be noted that it is not necessary that 
fi,m{f,p) is positive function, as long as those averages are obtained. In practice this means 
that initial conditions are randomly selected from \fi, m (r,p)\ and if N sets are chosen then 
the averages are approximately 

w ^ -2 ^ £f=i ifj x pjf sign [fi, m (rj,Pj)} 
<{rxp) >^ : 

E j= i sign [fi, m { r j,Pj)\ 

and similarly for the angular momentum. The function sign [fi,m(^j,Pj)] is the sign of 
fi >m (fj,Pj) for the set of initial conditions (fj,Pj). However, much more stringent condition 
is that the function fi tm (r,p) is stationary. The meaning of this can demonstrated on simple 
example. If the particle of mass M is free and the set of its initial conditions is (fj,pj) then 
its position after time t is r = fj + jit and velocity v = \%. From these values one would 
be able to calculate a new distribution function g(f,p), by requiring to have properties as 
previously described for the function fi tm (r,p). Stationarity requires that the two functions 
are identical for any time t. This restriction is very important, for obvious reasons, and it 
is sufficient to dismiss a large number of ad-hoc distributions. 

The most obvious starting point would be the wave function for a particle in a particular 
angular momentum state, which in the coordinate space is ij)i, m {r) = Xn,l( r )Yi,m(@i 0)j where 
Yi,m{Q, 0) is spherical harmonic. Likewise in momentum space the wave function is (fi, m {p) = 
u nt i(p)Yi tTn (9 p ,(j)p), where 9 P and <p p are angles of the momentum. One could then form 
the function fi >m (r,p) = \Xn,i{ r ) Y i,m{9, <P)\ 2 \uj nt i(p)Y l)m (9 p , 4> p )\ 2 from which initial position 
vector and momentum could be selected, except that neither the averages (D) are obtained 
nor the condition of stationarity for the function fi, m {r,p) is satisfied, which is easily proven. 
Therefore this simple idea fails, and this is because one basic difference between classical and 
quantum dynamics was not taken into account: classical dynamics is defined in the phase 
space, while quantum dynamics is either in the coordinate or momentum space. Simple 
product that was suggested for the distribution fi )m (r,p) only reflects this: it says that 
initial conditions in the phase space are determined as the product of the distribution from 
the coordinate and momentum space, as being independent. What is needed is formulation 
of quantum mechanics in the phase space, and then one would possibly be able to satisfy 
those two conditions. The earliest attempt in this direction was done by Wigner who 
looked for a function p(f, p, t ) with the property 

P(r,t) = \ij(r,t)\ 2 = Jd 3 pp(r,p,t) ; Q(p,t) = \^(p,t)\ 2 = J d 3 r p(f,p,t) (2) 

and obtained, what is known as, the Wigner function 

P(r,M = J^J d 3 qe 2 ^ h r(r + q,t)ij(f-q,t) (3) 

which is regarded as extension of quantum mechanics into the phase space. By defining 
it in this way inevitably results in non uniqueness of the extension H, because there are 
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a large class of phase space functions that are defined by the requirements (0). Another 
attempt to formulate quantum mechanics in the phase space is due to Moyal [|,|| whose 
main objective was to give it a sound statistical foundation. The main starting point is 
postulating observables as operators, and postulating that to each set of commutating ob- 
servables there is a set of non commuting. Together they form a complete set of observables. 
For these observable (operators) one forms a function from which the characteristic function 
of statistical theory for a given quantum state ip is obtained as a matrix element. By the 
Fourier transform Moyal obtains distribution which specifies to the Wigner function if the 
observables are coordinates and momenta. Common to both approaches is that the princi- 
ples of quantum mechanics are assumed, although Moyal analysis is more general. Can the 
Wigner function be used for the distribution /z, m (r, p)? It is not clear that the averages ([!]) 
are obtained, in fact as it will be shown later they are only partly obtained. Furthermore, 
it is also not clear that the stationarity condition is satisfied, however, the Wigner function 
is stationary under the quantum time evolution for the stationary wave function, which is 

gner function was used in various applications and its properties 
investigated [f7|-|2T|j, but its true significance is when connection with classical dynamics is 
sought. It is achieved by the standard assumption that classical dynamics is the limit of 
quantum when h — > 0. It can be shown in general (one dimensional problem is discussed for 



the moment) that the Wigner function satisfies the equation [19,20 



where 



and in this limit the inhomogeneous term vanishes. It also vanishes if the potential is 
quadratic, irrespective of this limit. The homogeneous equation is the Liouville equation, 
which determines time evolution of the phase space density p(r,p,t), and it is solved by 
classical equations of motion. 

Despite this, what appears, very important connection between quantum and classical 
dynamics, there are at least two remarks that one can make about this approach. It is 
desirable from understanding quantum-classical connection to learn more about foundations 
of quantum principles because with classical we are quite familiar. Strictly speaking that 
goal was not achieved because wether the Wigner or Moyal formulation of the phase space 
nothing can be learned about quantum principles because they are the starting point anyway. 
The second remark concerns the limit of the phase space density when H — > 0, which should 
also be taken and not only the limit for the inhomogeneous term. For the stationary states 
it can be shown that |13| 



p(r,p) =J(H-E )=S 



where Eq is a fixed energy. This form greatly restricts possible phase space densities from 
among those stationary solutions that are obtained from the Liouville equation. However, 
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it is consistent with the accepted view which is based on the correspondence principle: 
quantum goes over to classical solution for large quantum numbers (the proof of this is 
straightforward but not elaborated). 

Classical-quantum connection can be solved by starting from entirely classical principles 
p2| , |23|j , and as it will be shown allows much greater flexibility in application of classical 
dynamics to quantum problems. One starts from the Liouville equation in classical dynamics, 
and the argument on which one basis its use is quite straightforward: initial conditions 
for a particle are never accurately determined, in which case its precise trajectory has no 
meaning, only probabilities where it is found in the course of time. For example, to claim 
that position of Earth is accurately determined is nonsense, and therefore prediction of its 
position accurately in the next, say, million years is not possible, especially if perturbation 
from other planets is taken into account. However, it is possible to predict probability 
of finding it at certain position after (almost) any length of time. Therefore one starts by 
formulating classical mechanics as a statistical theory (it should be strongly emphasized that 
statistical does not imply many particles but probability for a single particle), in contrast 
to the attempt by Moyal whose principal aim was to formulate quantum mechanics as a 
statistical theory. In other words, one particle in the phase space is not a point but an 
extended density, which takes into account all uncertainties in determining its whereabouts. 
For a particle of mass m this equation is 

d t p(f, p, t) + - ■ V r p(f, p, t) + F ■ V p p(f, p, t) = (4) 
m 

where F is the force on the particle and p is its momentum. The index of the operator nabla 
designates the variable with respect to which the derivatives are taken. For a particular case 
when phase space density is stationary, which means that d t p(f,p,t) = 0, the solution is a 
function of the form p(r,p,t) = f(hi, /t2, hz, ■■■), where / is arbitrary function and hi are 
dynamic invariants of the classical equations of motion. One of them is Hamiltonian for the 
particle, and another are components of the angular momentum if the force F is centrally 
symmetric. There are other invariants but they will be mentioned later. These phase space 
densities, like any other solution of the Liouville equation, lack a very important ingredient 
in order to be regarded also quantum solutions. This ingredient is restriction on the possible 
phase space densities, and it is in the form of the uncertainty principle 

Ax Ap x y h/2 

for any Cartesian coordinate. The standard deviations Aa; and Ap are calculated from the 
phase space density. Mathematically speaking, selecting phase space densities according to 
this restriction is a very well defined problem and it is solved within the Fourier analysis. 
They are obtained in the form of the convolution 

p(?,p,t) = j d 3 qe 2i ^ n r(r + q,t)f(f-q,t) (5) 

which is also recognized as the Wigner function. The function / is arbitrary, but if the phase 
space density satisfies the Liouville equation then it satisfies the equation |24j 



indtf = ~Af + Vf (6) 
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where F = — W , and for potential it is not required to be harmonic (for a general potential 
parametrization of the phase space density is more elaborate). 

The common point between the Wigner and Moyal approach and the one that starts from 
classical principles is the function (|3]) or (|5]), however it is almost arbitrary (because there 
is no reason to choose it from many other functions) function in the former case but it is 
convolution in the latter. Apart from that common point there are fundamental differences 
between the two approaches, the former starts from quantum principles and the latter from 
classical. One manifestation of that is attitude towards the uncertainty relationship, which 
in the former case it is not considered a principle because it is derived from the other ones, 
while in the latter case it is assumed to be a fundamental principle from which the basic 
dynamics equation of quantum mechanics @ is derived. One could cite other fundamental 
differences between the two approaches, but one is that of simplicity. In the Wigner-Moyal 
approach one needs all postulates as in formulation of quantum mechanics, about whose 
number there is no general consensus, but in the approach from classical mechanics one 
needs only the postulates that defines it (i.e. there is no need for postulates other than the 
ones that one is already familiar width) and the uncertainty postulate. One could question 
the meaning of the uncertainty principle, but one could equally question the meaning of 
the observable-operator, or wave-particle dualism postulates in the formulation of quantum 
mechanics. Other advantages that one has by starting from the classical principles and 
not from the quantum phase space formulation of Wigner and Moyal, will be explicitly 
manifested in the following sections. In particular it is not clear how to derive the phase 
space density for the true rigid rotor, which is done in this work, by other than starting from 
the classical principles. 

There is one conceptual problem in the attempt to merge classical dynamics with the 
uncertainty principle. Without the latter solutions of the Liouville equation are the phase 
space probability densities, and as such they always have positive value. As soon as one 
imposes restriction in the form of the uncertainty principle this feature is lost and one often 
ends with the phase space density that has positive and negative values. This is the price 
it is paid by imposing that restriction, and its physics is justified on the grounds that the 
uncertainty principle makes it impossible to measure precisely the phase space probability 
density. One talks then about the phase space density, whose properties are exactly the 
same as for the probability density. However, all measurable quantities, e.g. probability 
density for the coordinates only, must have physically acceptable value. 

In the equation (||) one recognizes Schroedinger equation, the basic equation of quantum 
mechanics. Therefore the described steps in classical mechanics produce identical results as 
quantum mechanics, provided the initial condition for the phase space density is calculated 
from the parameterization (|5|). From these principles angular momentum for the three 
dimensional harmonic oscillator will be analyzed first. 

II. HARMONIC OSCILLATOR IN THREE DIMENSIONS 

Stationary phase space densities for a three dimensional harmonic oscillator will be an- 
alyzed if it is assumed that the average 

L= d 3 r d 3 p f x p p(f,p) (7) 
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has a given value. The starting point is the phase space density which is parametrized as 
where the functions if) are chosen in the form if) = r l Y t ^(9,(f)) R n ,i(r) = Yi t/Jl (x,y,z) R n j(r), 
where 9 and <f> are spherical angles and r is the radial coordinate. The indices n, I and p are 
integers (in further analysis the units are set in which m = h = 1 and the frequency of the 
oscillator is uj = 1), Y"; iM (#, 0) is spherical harmonic and 

R n , l {r)=N 1 F 1 (-n,l + ^; r 2 ) e"^ 

where N is normalization constant and ±Fi (a, b; z) is hypergeometric function. The phase 
space density is now 

pn,iA^P) = ^3 / ^qe^Y^ix + q x ,y + q v , z + q z ) R n ,i(\f+ q\) (8) 
Yi^{x -q x ,y- q y ,z- q z ) R Th i(\r - q\) 

which in general does not have simple explicit form, but it has few nice features. It is 
stationary, which means that its form does not changes with respect to the classical time 
evolution. In other words, if classical solution for the trajectory is 

r = r cos(t) + p sin(t) 

then 

p(r,p,t) = pnj,^ [f cos(t) -p sin(t),f sin(t) +p cos(t)] = p n ,i,p{r,p) 

This means that the phase space density is a function of the dynamic invariants of the 
harmonic oscillator. One set of these are elements of the energy tensor (| is omitted for 
simplicity) J5y = PiPj +XiXj ; i,j — 1,2,3 (the indices designate the Cartesian components 
x, y and z) , and the other are the components of the angular momentum 

; * 7^ j- 

Of course, any other combination of these basic invariants is possible, e.g. the total angular 
momentum squared. Which ones are present in the phase space density (§) are determined 
by explicit calculation, and the first few for the ground vibrational (n=0) state are given in 
Table I. The symbols represent the following quantities: E = p 2 + r 2 , E z = p\ + z 2 , L 2 = 
(iJPz — z Py) 2 + ( z Px — xPz) 2 + (xPy — UPxY and L z = xp y — yp x . These are dynamic invariants 
for the harmonic oscillator, and therefore the phase space density is indeed stationary. 

From the phase space densities one can calculate the total angular momentum ([?[) and 
its squared modulus from 

L 2 = J d 3 r d 3 p (fx p) 2 p(r,p) 

Their values are given in Table I and as expected angular momentum is the same as from 
quantum analysis: it has only the z component and its value is p. However, the angular 
momentum squared is not equal to 1(1 + 1), as expected from the quantum treatment, but 
differs by 3/2. The same is true for other than n = states, as shown in Table II. In fact 
the most surprising finding is that the states with I = 0, which are normally associated with 
the zero angular momentum, have the value 3/2 for the angular momentum squared. The 
question is where this discrepancy comes from? The simplest answer is that classical analysis 
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is not correct, because it surely must violate certain rules that are not consistent with the 
quantum mechanical ones. However, this answer is not correct, because if one writes the 
momentum squared as p 2 = p 2 +Pe +p$, where the components of the momentum are given 
with respect to the vector f, then its average for the function tp = 0(0, <j))R(r) is 



< 



1 r 

P 2 >= ^3 J d ^ r d Pr d P0 d P<P (Pi +Pd+ Pi) 

J dq r dq e d q<t> e 2 <^+^+*w)e*(0 + , + )i2(r + )6(0_, <j>-)R{r. 



where 



r± = J(r ± q r ) 2 + q$ + q\ ; cos 6± 



z-(f±q) (r ± q r ) cos6±qe sini 



\r±q\ yj(r±q r ) 2 + q 2 e+ql 



(9) 



?i<P± 



(f± q) ■ (x + iy) _ (r ± q r ) sin 9 =F qe cos 9 ± iq^ H 
\f ± q \ sin 9± 



\f ± q\ sin ± 

After a straightforward, but lengthy, calculation of the integrals one obtains 

,d\ 1 d ( . _ d \ 1 9 2 



<p 2 >=- J d 3 r i>*(r)^ 



d 
dr 



+ 



dr J sin(9) 06 



sin{9) 



86 



+ 



sin 2 (6) 



which is the correct answer for the kinetic energy operator (up to a pre-factor, which was not 
taken into account). In the standard interpretation the angular part is then associated with 
the angular momentum squared operator, which indeed gives zero for angular momentum 
squared for the I = states. Therefore the procedure of deriving the kinetic energy operator 
is correct, but then the question is where the analysis that produced results in Tables I and 
II is inconsistent with the quantum interpretation? To answer this question one explicitly 
calculates the average of p^ = p 2 d +p^, which is the classical angular part of the momentum 
squared. For simplicity its average will be calculated for / = state. 
By definition 



< 



1 r 

Pn>=j^J d * r d Pr d Pe d P<p {fe + pf) 

dq r dq e d Q<p e 2 < p ^ +pM ^) R(^r + q r ) 2 + q 2 + q 2 )R(y/(r - q r ) 2 + q 2 + q 2 ) 



where integration in p r and q r is calculated first, then p 2 d and p 2 are replaced by derivatives 
in the variables qg and q^ that act on the exponential function, respectively, and after partial 
integration in the same variables one obtains 



<Ph >-- 



16n 3 



d 3 r dpg d P<t> f dq e dq^ e 2 < Peqe+p ^ (d 2 qe + & t 



I r 2 _l „2 _l „2 



r + Qe + 



After integration in the variables pe and p<f> the only non zero contribution is 

< pi >= -— f d 3 r R{r) -d r R{r) 
Air J r 

and it is not zero. In fact the average of the angular momentum squared (which is L 2 = r 2 p^) 
is 
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< L 2 >= — / dr r 3 R(r) d r R(r) = - 
Jo 2 



which is precisely the value in Tables I and II. Furthermore the result is independent of 
the radial function, which indicates that this is a universal number. In order to check the 
correctness of the result one calculates the average of the radial part of the momentum 
squared, and the result is 



d 2 R{r) + -d r R{r) 
r 



< Pr >= ~ J d 3 r R(r) 
Together with the angular part one obtains the radial part of the kinetic energy operator 

T r = -^d r (r 2 d r ) 

which is the correct answer. Therefore, there is no discrepancy between quantum expression 
for the kinetic energy and classical treatment in this work, but the difference is in the 
interpretation of its fragments. According to results in this work one should interpret the 
radial kinetic energy operator as 

t = -dl - -a 



r 



while 



1 d ( d \ 1 & 

sm(9)— + 



sin(6)d6 \ v ' 89 J sin 2 (6)d<p 



is operator for the angular momentum squared. 

Further confirmation of this comes from calculating the average angular momentum 
explicitly. For I = states it is given by 

< L >=< L z > z = J d 3 r dp r dp e d P<p r P4> sin(^) J dq r dq e d q<j> e 2 < p ''^+ p ^ s +^^) J R(r + ) J R(r_) 

which is easily shown to be zero. This is no contradiction with the previous finding, because 
zero of the angular momentum is result of cancellations of the contributions from two signs 
of rather than the modulus of the angular momentum being zero. One can therefore 
calculate angular momentum by averaging over only one sign of p^, because the average 
over the other gives the same result but with the opposite sign, and the sum-total is zero. 
This average is 

< L+ >= J d 3 r dp r dp sin(0) d P4> r V<i> J dq r dq e d q<t> e 2i ( Mr+Peg9+p ^) R(r+)R(r-) 
and after evaluating most of integrals one gets 

By using the relationship 
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PCX) 

Jo ^ 



(2?*, 



where designates the principal value of the integral, it is obtained 

< Lf >= J #r rsin(^) % J dq, R{^r 2 + + qj) -j== = 

The principal value can be omitted because the integrand is not singular, in which case 

< L+ >=- [dr r 2 R 2 (r) = - 
2 AJ 4 

The value of the z component of the angular momentum, which gets contribution from the 
space > 0, is independent of the radial function R(r), the same result as before. 



III. RIGID ROTOR 



Rigid rotor in classical physics is a very well defined object, but one could argue that 
it is also in quantum. Unfortunately, as it will be shown, the two definitions have not the 
same meaning. Intuition, which is essentially classical, defines rigid object as the one which 
does not have breathing (radial) motion. Objects that appear to qualify as the rigid are 
diatomic molecules (at least many of them), which is supported by the approximations in 
their quantum description. Schroedinger equation for the nuclear motion in the diatomic 
molecule is 



%Ti~z — = — - — 
dt 2m 



1 d 



+ 



1 



d 



r 2 dr \ dr J r 2 sin(6) d6 



sin(9) 



d^\ 1 d 2 i> 

~d9J + r 2 sin 2 (9) d<j) 2 



+ V(r) i) 



where m is the reduced mass of the diatomic molecule and V(r) is the internuclear potential. 
Typically this potential has a deep minimum displaced from the origin by the distance r (the 
bond length), which for the modelling purpose is approximated by |w 2 (r — r ) 2 . Ground 
state wave function in this potential has the width a, and its relationship to the bond length 
is r ^> a. Under this circumstance the radial coordinate in the angular part of the kinetic 
energy operator is approximated by the fixed value r , in which case the wave function 
factors as tp — ^R(r)Q(9, (f),t)e~ lEot / h , where E is the ground state energy of the molecule. 
The angular function then satisfies the equation 



ih 



dQ(9,(f),t) 
~di 



T rot e(6,<i>,t) + w(e,<i>,t)e(e,<i>,t) 



(10) 



where 



T-rot 



h 2 



2mrQ 



sin{9)d9 ( Sm{9) d9 / 



+ 



d 2 



sin 2 (9) 



and W(9, 0, t) represents interaction of the molecule that depends on the angles only (say a 
dipole in the electric field). 
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The unique feature of the equation (|T0|) is that it depends only on the angle variables, 
orientations angles of the molecule, and in this sense it is called the equation for the rigid 
rotor. The operator T rot then represents the kinetic energy operator for the rigid rotor, 
which has discrete spectrum 



h 2 

"o 



T r ot Y ltm {e,<j>) = - — ji(* + i)*WM) ; / = 0,1,2,.. m = -/,-/ + 1,....,/ 

ZTYlVn 



where Yz, m (0, 0) are spherical harmonics, the eigenfunctions of the rigid rotor. However, this 
is far from being true, if one has the image of the classical rigid rotor in mind. Although 
in the previous derivation dynamics in the angle variables is de-coupled from the dynamics 
in the radial this is not sufficient to call the molecule a rigid rotor. It is also necessary to 
prove that the radial component of the momentum for the relative motion of the two atoms 
is de-coupled from its angular ones, and to show this one needs to calculate the momentum 
space wave function. This function for a stationary rotational state of the "rigid rotor" is 

d\ -R(r)Y ljm (6,4>) e ipV = (2irf 2 i% m (9 p ,<f> P )^ / dr r 1 ' 2 R(r)J l+1/2 {pr) 
r y/p Jo l^) 

where the expansion of the plane wave in the spherical coordinates was used. De-coupling 
of the radial component p r from its angular ones pe and p^ is therefore not possible because 
P — \/Pr + Pe + P%i while the spherical angles 9 P and 4> p are related to the same components 



in a complicated way. Therefore it cannot be assumed that the theoretical model that was 
just described represents rigid rotor. 

Analysis in the phase space is more revealing, but it is somewhat approximate. It was 
mentioned that around r = r the potential is very well approximated by harmonic, in which 
case the entire wave function for the diatomic molecule is 

where 1/r coefficient was assumed to be constant. The phase space density is then 

N 2 I" o r,-^ (\r + q\-rn) 2 (|r-g|-r ) 2 

Pl,m(r,P) = ^/Ae*"^,^)^-,^ 5^ sr*- 

where the angles 9± and <p± are defined in The rigid rotor model assumption implies 
that q << r, and if spherical coordinates are used for the vector q = q r f + q^9 + q^cj) then 
the exponent of the harmonic oscillator function is approximately 

(\r + q)-r ) 2 (]r- q\ - r ) 2 ^ (r - r ) 2 $ + {«0 + 4f_ 6 

la 2 la 2 a 2 a 2 \ r\ J Aa 2 rl 1 ' 

where in the expansion coefficients, except in the leading one, it was set r = r$ . By 
recalling that this expansion is in the exponent it follows that the range of \q r \ within which 
it significantly contributes to the phase space density is of the order a, however, the range 



of Jqg + q 2 is of the order ^ar . This means that \q r \ << y q 2 + q 2 in which case the phase 
space density is approximately 
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2r a 2 N 2 _(i^ol!_„2 a 2 



IT 

roc 



5/2 



f° dq e r dq^ e *^(™+r«*)-(ti+#) 9, 0, 0, V^qe, V^q^) 

J— oo J—oo 

where the angular function G/ m is shorthand for the product of two functions Yi m (9±, <j)±). 
Dependence of the angular function on the radial variable q r is also neglected, in which case 
the phase space density parametrizes as a product of the function for the radial variables 
and the function for the angular. The angular phase space density is now defined as 

f%M<l>,Po,P+) = ^ Idqe dq+ e W2^(i* W ^)-(#H$) a e ^(r , 0, <f>, 0, y/2^q 9 , y/2^q^) 
71 J (12) 

and appears to describe the rigid rotor because it is independent of the radial variables. 
However, this conclusion is false because the phase space density depends on the radial 
parameter a, and in this respect it depends on the radial dynamics. Because of this feature 
the rigid rotor will be called "soft". 

Average of the angular momentum squared has now the value 

<l 2 >= J (%JP,<I>,v9,pM {pI+pI) = Ki + i) 

where the constant term 3/2 that was obtained in the previous section is missing. This is 
expected because coupling with the radial motion was neglected. Another property of the 
angular phase space density is 

/oo roc 
dpe / d V4> plJBA.P^) = 6 i)m (r ,M, 0,0,0) = \Yi, m (0 , <j>)\ (13) 
-oo J—oo 

which shows that it correctly describes the essentials of the angular momentum. Its explicit 
form for I = m = is given by 

Po,o(0> <f>,Pe,P<p) = -^rJ Q d 1 Q J o (2V2ar pg) e"" 4 

which does not have explicit analytic expression, but its shape is simple, as shown in Figure 
1. For the rotational state I = 1 and m = the angular phase space density is 

Ar 2- P e ' 4r 2"^)/ dqqJ (2V2a7'oPq)e- qi 



( 2 „ 1 n2 C °S 2 8 

~ cos e + —,d 2 + 



where e = a/r , and in the last step only the terms up to the order e where retained. 
Similarly the phase space density for I — 1 and m — 1 is 

(sin 9 1 1 sin 2 9 \ f°° / \ 

Sin2 9 + 1^ dp * + ~&f ve + ~&f V4 > + Jo d<1 qJ ° ( 2v/ ^ W ) 6 
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These phase space densities, except Pq , are time dependent, which means that they are 
not functions of only the dynamic invariants of the rigid rotor, e.g. p 2 = p 2 e + p 2 and p^ sin^. 
Thus for example the terms in Pi {9, 4>,pe,Pcf>) that cannot be represented by invariants are 



Pxo{Q,<t>,Pe,P<j>) 



roo , . 4 COS 2 9 I QjTq Z" 00 / \ 4 

cos 2 9 J dq qJo (2\/2aropq) e~ q ~\~2~ J ^ [pV^/ropqJ e~ q 



Explicit time dependence of the phase space density is obtained from the time dependence 
of the angle 9 

a a tpo p° e . . tpo 
cos 9 = cos 6*0 cos sm 9 sm 



mro po mro 

where 9q, Pg and po are initial values of these variables. In the phase space density cos# is 
replaced by 

a a tP , Pe ■ a ■ tp 
cos 9 t = cos 9 cos 1 sm 9 sm 



mro P 

and the time dependence of the probability density (|13|) is calculated from 

dpe \ dp^ p 1>0 (9,(j),Pe,Pi) ~ / dp g dp<p cos 2 9 dq qJ [2y/2ar pq) e 

-oo J —oo J —oo J —oo JO 

= 1 + cos 9 . 3 cos 9 — 1 )d t I dq 

32ar V J 8^2^ V J J H 



</ 4 



2am 2 " 



The first part is time independent, while the second goes to zero after the time interval 
t > mr ^2ar , which is typically of the order 10 -13 — 1CT 14 sec for diatomic molecules. 
Therefore the probability density starts as flllf ) but its limiting value is constant, but not in 
the form of the squared modulus of the spherical harmonic. 



IV. TRUE RIGID ROTOR 

Previous discussion revealed great difficulty in formulating the concept of the rigid rotor 
in quantum theory. This was manifested as inability to formulate the phase space density 
that involves only the parameters for the rotational degrees of freedom. It does not help to 
take the limit a — > in ( |T2] ) because that would imply infinite dispersion of the variables 
pe and p,p, which only reflects the fact that the radial and the angular components of the 
momentum are interrelated. There is, however, a way of formulating the true rigid rotor, but 
should be done by following formulation of quantum mechanics as suggested in Introduction. 
One starts from the Liouville equation in the spherical coordinates, which for a free particle 
is 

a PO A , P<t> A , Pr a 

dtp o e p H -r—^o^p H d r p - 

mr mr sm9 m 

1 (p^cosO \ 
— ■ a +PrPe o Pe p 
mr \ sm.9 I 
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and the rigid rotor assumption implies that the phase space density is r and p r independent. 
This means that the Liouville equation for the rigid rotor is 



dtp d e p H '— -d^ + 

mr mr sm9 



cos 9 



mr sin 9 



d H p - P4, d Pe p 







(14) 



where r is constant. It can be easily verified that if the phase space density is a function of 
the form p(9,(j),pg,p ( f > ,t) = F{p 2 e + p^p^ sin 9), where F is arbitrary function, then it is time 
independent. Additional requirement is that the phase space density should be in accordance 
with the uncertainty principle, which is achieved by straightforward generalization of the 
rule that was used before. The phase space density is therefore parametrized as 

= -J d 2 qe 2 < p ^ +p ^) e*(9 + ,<p + ) 0(0_,0_) (15) 

7T Z J 

where the angular functions will be determined for a particular case when the solutions 
of the Liouville equation (14]) are stationary, i.e. d t p = 0. The relevant variables were 
defined in (g). Parametrization (|T5| ) is replaced in the Liouville equation ([14]), and by using 
transformations of the kind 

p e p(9,<p, Pe ,pt) = -Ljd 2 qe 2 < p ^ +P ^)d qs [e(0_,^_) e*(# + ,0 + )] 

p 9 d PeP (0A,P8,Pt) = - / d 2 qe 2 < p ^ +p ^)d qe [q e 0(0_,0_) 0*(0 + ,0 + )] 



one obtains, after lengthy simplifications, that the stationary solutions of the Liouville equa- 
tion (0) satisfy 



where 
fi= & 
and 



^/ dV ^ +w )_±_( /l + I /2 ). 

in H V / sin fl v ~ sin H V 



h = _ 7r( cos ^+ + cos6>_ 
2r 



(16) 



1 fl2 

sin 2 ^ ^ 



sin 6L 



sin 6* i 



COSt^i — costM — n — 

2 V + ' \ sinfl„ sinfl. 



, sin 9 



e*fy_e 0^0* 



sin' 9. 



+ 



sin' 0, 



where q 2 = q 2 + q\- It is implied that is a function of the variables 9- and <p- while 0* is 
a function of 9 + and 0+. 

If the function /2 is neglected for the moment then the condition (|1(J) implies that the 
angular function satisfies the differential equation 



1 

sin# 



d e {sin9d e Q) 



sin 2 9 



die = A 
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where A is a real constant. In this equation one recognizes the equation for the spherical 
harmonics, where A = —1(1 + 1). However, the function f 2 cannot be neglected, and therefore 
the angular functions are only approximately the spherical harmonics. However, a very useful 
feature of the function f 2 is that in the limits qe — > and q^ — > it is equal to zero, in 
which case the angular functions are exactly the spherical harmonics. This means that when 
the phase space density is integrated over the momentum variables the resulting probability 
density should be the squared modulus of the spherical harmonics, i.e. 

P(M) = Jdpe Jdp^p(e,(t>, P6 ,p^) = |6(M)| 2 = \Yi, m (d,<i>)\ 2 (17) 

The choice of the spherical harmonics for the angular functions is that the phase space density 
is approximate, which is manifested as being time dependent, i.e. it is not a function of only 
the dynamic invariants of the rigid rotor. This is the price that is paid by neglecting the 
function f 2 in the equation (|I7J). Inclusion of this function results in the phase space density 
that is a function of only these invariants, and this fact is used as the procedure to find a 
proper phase space density. This is best demonstrated on one example. One particular case, 
however, has exact solution, and this is when the angular function is constant. In this case 

PD,o(M,P»,P*) = i/Ae 2i ( M8+w ^ = ^W%) 

The example that will be analyzed in more details is when the angular function is the 
spherical harmonic Yi^(9,<fi). The phase space density (p~5|) is 

p(9,cj ) ,pe,p <t> ,t) = -1 [ d 2 qe 2 <™° +p ^) Y? (9+,<p+) Y 1>0 (9-, <j>J) 
and its explicit form is (the constant r is fixed to unity) 



4jr 3 J * 1 + q* 

where K n (x) is modified Bessel function of the second kind. This phase space density is not 
stationary, because the term 



tj Dill , 

Pt = ~^T K i^P) 



3 sin 2 9 

ApTT 2 

is not a combination of the dynamic invariants. However, if p^ is replaced by p cos a then 



2W da (l-2^) = 
o \ p l I 



which means that these two terms in the phase space density can be omitted without in any 



way modifying the integral (|I7), the value of the total angular momentum squared 



< L 2 >= J dn J d 2 p p 2 p(9,(fi,p e , 
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and the angular momentum (its z-th component) 

< L >= z J dQ J d 2 pp^sia9 p(9,(f),pg,p ( j ) ) =< L z > z 

Therefore 

3 ( jy^ sin^ \ 3 
(h,o(0,<l>,Po,P*) = ( 1 - ^2 — J J ^o(2p) - — 5(p e )5(p (l> )sm 2 e 

is time independent and represents phase space density for the rigid rotor in the state with 
the angular momentum squared < L 2 >= 2 and angular momentum < L z >= 0, while the 
probability density ( |TTD is P(8,<p) = \Yi t o(6, 4>) f- Similarly the phase space density for 
I = 1 and m = 1 is 

MM,*,*) = ^ (l + K (2p) + 3 J^ Kl (2p) + ±- (-2 + sin 2 *) 5(p e )5(p,) 

with the property that P(6, <p) = \Y ltl (9, </>)| 2 , < L 2 >= 2 and < L z >= 1. 

The phase space densities for the angular momentum states (2,m) ; m = 0,1,2 were 
calculated as additional example. They are associated with the angular momentum squared 
value < L 2 >= 6 and the angular momentum < L z >= 0, 1, 2, respectively. They are 

15 / plsin 4 #\ , s 45p / p 2 6 sm 2 6\ 2 , x 

(8 - 24 sin 2 + 27 sin 4 fl) Sfatffo) 

15 / p|sin 2 fl _ , 2pUin 4 9\ , s 
MM,*,*) = ( - 1 + 2*sin0 + ^ 2p 3 sin 3 # - ^ j K (2p) + 

15P ^ _ 2 ^ _ Pl^l) Ki{2p) + W / in2 , _ 3 gm4 A 



4tt 2 \ p 4 p A J y r/ 16tt V 4 

15 / »;isin 2 # 7rj,sin 3 # pi sin 4 0\ 



15p ^ _ + 6 g|sin^ + 4p3 sin 3 S + g^sin^ y + 



167T 2 y p 2 p 2 3p 4 p A 

^L(l-^ 2 6 + l^6) S(pe)S( P4> ) 



In this way phase space for the true rigid rotor was defined. It is "true" because only 
the variables that are relevant for such object were introduced. 



V. DISCUSSION 

Analysis of angular momentum in the phase space was made, and in this context rigid 
rotor was discussed. Perhaps one of the most intriguing finding is interpretation for the 
partition of the kinetic energy operator, which is in considerable disagreement with the 
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standard one. Part of what is considered to be radial kinetic energy is in fact contribution 
from the angular momentum operator, despite the fact that it only contains radial vari- 
able. This finding does not come as a surprise if one makes the following observation in 
the traditional classical mechanics. Given spherically symmetric probability distribution of 
coordinates P(r) and momenta Q(p) for a free particle (say it was obtained my measuring 
its position and momentum) the average kinetic energy is not zero despite the fact that the 
momentum is. Likewise, the average angular momentum is zero but the average modulus 
squared of it is not because 

< L 2 >= J d 3 r J d 3 p r 2 (p 2 d +pl)P(r)Q(p) ^ 

Introducing the uncertainty principle does not changes this fact, except that instead of this 
average having arbitrary value it has a fixed and equal to 3/2. Therefore a particle has always 
non zero modulus of angular momentum, and its minimal value is fixed and independent 
of the phase space density. It is like saying that particle always carries a minimal intrinsic 
angular momentum, but in the way it is described it never manifests itself. For all practical 
purpose this finding is immaterial, because standard interpretation is self sufficient, but it 
becomes evident when classical modelling is attempted. 

From the phase space density for the rigid rotor one expects to obtain momentum space 
probability density, which is given by 

Q(P) = J d3r P(r,p) 

however one should be careful about components of the momentum variable. Throughout 
the paper the components with respect to the vector f were used, because they are natural 
when angular momentum is analyzed. This is because pe and p^ are components along the 
appropriate angular unit vectors that are perpendicular to the vector r and hence directly 
proportional to the angular momentum (they can be called radius vector components of the 
momentum). However, in the momentum space one works only with the components of the 
momentum vector, which for the mentioned components this is not the case. If p is modulus 
of the vector p while P and <ft p are its spherical angles then 

p r = P (cos cos 6 P + cos [4> — <p p ) sin sin P ) (18) 

Pe — P ( srn cos $p — cos [0 — 4>p] cos ^ s i n ®p) 
p^ = -p sin [<p - 4> p ] sin p 

which explicitly shows that the radius vector components of momentum are a mixture of 
momentum and radial vector spherical coordinates. Therefore in the phase space density 
they must be replaced by fll8|) and then the integration over the spatial coordinates per- 
formed. Indeed for the "soft" rigid rotor it can be shown that one obtains for the momentum 
distribution the square modulus of (PD, but it is not clear what the outcome would be for 
the true rigid rotor. For the latter it is required that p r = and yet the expression for the 
phase space density would be a function of p and P without an obvious restriction of that 
kind. However, explicit expression for the momentum probability is not of importance, it 
is important to be able to calculate the averages. For example the average of the square of 
the Cartesian component p x , say with the phase space density p 2 ,o is given by 
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22tt 



/ ZZ' 
dVL dp e dp^ [p^ sin <p + p e cos 6 cos <p} 2 p 2 ,o (0,<j>,Pe,P<i>) = — 

where p r = was set. 
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FIGURES 

FIG. 1. Phase space density for "soft" rigid rotor, for the angular momentum indices l=m=0 
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TABLES 



/, n 


Po,l,n ■ K 3 e E 


<L> 


<L 2 > 


0,0 


1 





0+| 


1,0 


-1 + 2E Z 





1*2 + | 


1,1 


—l + E — E z + 2L Z 


1 z 


1*2 + | 


2,0 


1 - | (E + 3E Z ) + l(E- 3E z f + §L 2 - 4L 2 





2*3 + | 


2,1 


(1-2E Z )(1-E + E Z -2L Z ) 


1 z 


2*3 + | 


2,2 


-l + \[E-E z - 2(1 - L z )f 


2 z 


2*3 + | 



TABLE I. Phase space density for the ground vibrational state of harmonic oscillator for the 
first few values of angular momentum numbers. Appropriate angular momentum and its modulus 
squared are given. Definition of the variables is given in the text. 





Pa,i,n ■ 7r 3 e B 


<L> 


<L 2 > 


0,0 


1-| (4E - 2E 2 + 8L 2 ) 





0+| 


1,0 


-1 + t 


( 2E-E 2 + 9E Z - 8EE Z + \ 
y 2E 2 E Z + 12L 2 - 8E Z L 2 - 8L 2 Z ) 





1*2 + | 


LI 


-1 + 1 


1 13E -WE 2 + 2E 3 -9E Z + 8EE Z - \ 

2E 2 E Z + 16L 2 - 8EL 2 + 8E Z L 2 + 10L Z - 
v 8EL Z + AE 2 L Z - IQL 2 L Z + 8L 2 J 




1 z 


1*2 + | 



TABLE II. Phase space density for the first excited vibrational state of harmonic oscillator 
for the first few values of angular momentum numbers. Appropriate angular momentum and its 
modulus squared are given. Definition of the variables is given in the text. 
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